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Abstract — In this letter, we propose a new molecular modula- 
tion scheme for nanonetworks. To evaluate the scheme we intro- 
duce a more realistic system model for molecule dissemination 
and propagation processes based on the Poisson distribution. We 
derive the probability of error of our proposed scheme as well as 
the previously introduced schemes, including concentration and 
molecular shift keying modulations by taking into account the 
error propagation effect of previously decoded symbols. Since 
in our scheme the decoding of the current symbol does not 
depend on the previously transmitted and decoded symbols, we 
do not encounter error propagation; and so as our numerical 
results indicate, the proposed scheme outperforms the previously 
introduced schemes. We then introduce a general molecular 
communication system and use information theoretic tools to 
derive fundamental limits on its probability of error. 

Index Terms — Nanonetworking, Molecular Communication, 
Modulation, Poisson Distribution. 

I. Introduction 

NANONETWORKING promises new solutions for many 
applications in the biomedical and industrial fields. This 
new paradigm utilizes various methods to communicate in- 
formation between nano-scale machines [1|. A promising 
communication method at this scale, is molecular communica- 
tion. Different molecular communication schemes have been 
proposed that can be categorized based on their effective range 
of communication, namely short, medium or large 0. 

In this paper, we focus on a diffusion based communication 
system for short and medium range nanonetworks. In this 
system, molecule concentration represents information. The 
transmitter emits molecules whose type and intensity depend 
on the current input symbol. These molecules propagate 
through the environment, with a portion hitting the receiver's 
surface. Receptors located on receiver's surface form chemical 
bonds with the incoming molecules, initiating a process that 
eventually results in decoding of the transmitted information. 

Different channel models for molecular communication 
have been developed and their related channel capacities 
have been evaluated in (3) -(6). Authors in [7 | have proposed 
two modulation methods called Concentration Shift Keying 
(CSK) and Molecular Shift Keying (MOSK). In these schemes, 
information is encoded in molecule diffusion rate and molecule 
type, respectively. Both CSK and MOSK suffer from interfer- 
ence caused by molecules from previous transmissions, hitting 
the receiver after a long delay. To control this interference in 
CSK (where only one molecule type is used), the decoding 
of the current symbol is adapted to the last decoded symbol. 
Although this helps to manage the interference, it results in 



an error leakage to the current symbol if the last symbol is 
decoded incorrectly. Using multiple molecule types, MOSK is 
more resistant to interference than CSK but it requires complex 
molecular mechanisms at both the transmitter and the receiver 
for message synthesis and decoding Q. 

In this paper we propose a new modulation method based 
on the idea of using distinct types of molecules for consec- 
utive slots at the transmitter, thus effectively suppressing the 
interference. Although simple in its formulation, the method 
outperforms the existing ones. To evaluate the schemes, we 
propose a new system model for the molecule propagation 
process and the channel based on the Poisson distribution. 
This model is more realistic and amenable to analysis than 
the one proposed in [7| based on the binomial distribution. 
We evaluate and formulate the probability of error for our 
proposed scheme as well as the previous ones and compare 
the results. Specially, we take into account the dependency 
of decoding of the current symbol to the last transmitted 
and decoded symbols to compute these metrics. Numerical 
results indicate our scheme has a lower probability of error 
compared to CSK and MOSK. To understand the efficiency 
of our proposed scheme, we use information theoretic tools to 
derive a fundamental limit on the minimum probability of error 
of an arbitrary scheme. A comparison indicates that although 
our scheme outperforms the existing ones, its error exponent 
(in terms of transmission power) is larger than our fundamental 
lower bound, suggesting room for further improvement. 

The rest of the paper is organized as follows. Proposed mod- 
ulation scheme and system model are described in Section II. 
In Section III performance of different modulation techniques 
are discussed. Also in this section a more general molecular 
communication system is introduced for which a lower bound 
on probability of error is computed. Lastly numerical results 
are presented in Section IV. 

II. Proposed Modulation Scheme and System 
Model 

A. Previously introduced modulation schemes 

The Concentration Shift Keying (CSK) scheme is inspired 
by the Amplitude Shift Keying (ASK) used in the classical 
communication. In this scheme, symbols are encoded in the 
number of molecules that the transmitter diffuses; this is called 
diffusion rate or transmission power. For example, to represent 
b bits, 2 b different molecule diffusion rates are utilized. At the 
receiver, the concentration (number of the received molecules) 
is compared with decision thresholds in order to decode the 
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transmitted symbol. The molecules from the previous symbol 
in the last time slot may hit the receiver in the current time slot. 
Thus, to reduce the effect of the interference in this method, 
decision thresholds for the current symbol need to be adapted 
to the last decoded symbol. 

Molecular Shift Keying (MOSK) scheme resembles the 
orthogonal modulation in the classical communications. To 
transmit b bits per time slot, 2 b different types of molecule 
are utilized. The transmitter releases a specified number of 
molecules of the type corresponding to the current input 
symbol. To decode the transmitted signal at a certain time 
slot, the receiver looks for a unique molecule type whose 
concentration exceeds a certain threshold r. An error occurs 
if the concentration of none of the molecule types or more 
than one molecule type exceeds r. Although interference due 
to the previous transmission in MOSK is less than that in the 
CSK, but it still exists. Also using different molecule types 
complicates the hardware structure at both the transmitter 
and the receiver as the number of molecule types increases 
exponentially with b. 

B. Proposed modulation scheme 

The performance of both the MOSK and CSK schemes is 
degraded by the interference caused by molecules of previous 
transmissions. In addition CSK suffers from propagation of 
errors since the thresholds for each time slot is adapted to the 
previously decoded symbols. Therefore, if the previous symbol 
is falsely decoded, its error leaks to the current symbol. 

Our scheme borrows ideas from CSK and MOSK. Having 
molecules of types A% and A2, the transmitter uses type Ai 
in odd time slots, and type A2 in even time slots. The use 
of two molecule types resembles MOSK but the difference is 
that the molecule type is not used for signalling. To convey 
information in each time slot, different diffusion rates are 
used (similar to CSK). Thus, to match each symbol to b 
bits, 2 b different propagation rates are utilized in each time 
slot. As the molecule types are different in two subsequent 
time slots, (i) the previous symbol interference has almost 
disappeared (ii) the decision threshold of the current symbol 
is independent of the last decoded symbol. Since the data is 
not encoded in the molecule types (but in the concentrations), 
the number of molecule types (and as a result the complexity), 
does not increase with b. We refer to this scheme as Molecular- 
Concentration Shift Keying (MOCSK). 

C. A Poisson model for molecular communication 

In this section we develop a model for a general molecular 
transmission. Later we use this model to compare probabilities 
of error of MOSK, CSK and MOCSK. It is assumed that 
the symbols are transmitted in equal time slots called symbol 
duration, t s . Generally the signal is decoded based on the 
concentration of molecule types received during a time slot. 
We assume that in a time slot three sources contribute to the 
received molecules of a specific type: (i) molecules due to 
transmission in the current time slot (ii) the residue molecules 
due to transmission in the last time slot (the residue molecules 
from two or more earlier time slots are ignored) (iii) the 



molecules from other sources as the noise of environment. 
More precisely, each molecule either misses hitting the re- 
ceiver, or hits the receiver in the current or next time slot. Let 
the hitting probabilities in the current and next time slot be P\ 
and P2, respectively (P1+P2 < 1). We assume intra-molecule 
collisions have little effect on the molecule's movement. 

For each molecule type, the transmitter has a storage. Every 
storage has an outlet whose size controls how many molecules 
can exit the storage to be diffused in the environment. The 
distribution of the number of molecules that exit from the 
storage in a time slot is Poisson. This follows from the fact 
that the number of molecules in the storage is large and 
probability of each passing through the outlet is small. We use 
Poisson(X) to denote the number of molecules exiting the 
storage (and diffused in the environment). The parameter X is 
called the diffusion rate and is determined by the size of the 
opening of the outlet. Parameter X itself is a random variable 
that is decided by the transmitter based on the message it 
wishes to transmit. 

Assume that the transmitter has chosen diffusion rates of 
X\ and X 2 of a molecule type like A in two consecutive time 
slots in a modulation scheme. Also, consider the environmet 
noise independent of the source signal as a Poisson noise 
with parameter Ao- Considering hitting probabilities P\ and 
P2, and using the thinning property of Poisson and the fact 
that sum of independent Poisson rv's is itself a Poission, the 
number of type A molecules received within the current time 
slot (denoted by Y^ A ') follows a Poisson distribution with 
parameter A where A = P\X\ + P2X2 + Ao. This Poisson 
model has advantages over the binomial model considered in 
0: (i) The Poisson model is more realistic and practical in 
comparison with the binomial model. For instance, in binomial 
model the transmitter needs to count the number of molecules 
that it diffuses in comparison with the Poisson model in which 
the size of the storage's outlet is controlled (without counting 
the exact number of molecule passing through the outlet), (ii) 
The Poisson model is easier to work with analytically because 
of the nice properties of Poisson (e.g. thinning property, etc.). 

D. Decoding probabilities for CSK, MOSK and MOCSK 

In the CSK method, only one type of molecule for instance 
A is used in all time slots. So, the decoding is made by 
comparing the number of received molecules, i.e. Y^ A \ with 
the thresholds determined based on the last decoded symbol. 

In MOSK, as information is encoded in molecule type, if 
two subsequent symbols are the same, only one molecule type 
for instance Ai is received, otherwise two molecule types for 
instance A\ from the current symbol and A 2 type from the 
previous symbol are received. In the former case, the symbol 
is decoded correctly if the total sum of received molecules, 
i.e. Y( Al \ exceeds the threshold r. In the latter, decision is 
made based on Y^ Al \ y(^ 2 ) and the threshold r. 

In the proposed MOCSK scheme, as stated above, 2 differ- 
ent types of molecules in 2 subsequent time slots are used. 
As the receiver knows the type of the current symbol, it needs 
only to compare the number of the molecules of this type with 
the thresholds to decode the current symbol. We now compute 
the average probability of error for CSK, MOSK and MOCSK. 
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As in 0, we only consider the class of simple symbol-by- 
symbol decoders that are practically appealing. Comparisons 
are made consistently for decoders in this class with the same 
transmission powers. 

For binary CSK (BCSK) using molecules of type A, if the 
current input symbol is s c = and the last transmitted and 
decoded symbols are s p and ~s p respectively, the probability of 
successful decoding is computed as 

P s e \s c ,s p ,S p & = 0|sc = 0, Sp ,s p ) = P{Y {A) < r 5p ), 

where s c is the decoded current symbol and 73- is the decision 
threshold corresponding to the value of s p . Note that in this 
equation, s p is not necessarily equal to the previous transmitted 
symbol, i.e. s p . In fact, the decision threshold for the current 
symbol is determined based on the previous decoded symbol 
regardless of being decoded correctly or not. 

For binary MOSK (BMOSK), assume molecules type A\ 
and A 2 are sent for and 1, respectively. For the current 
transmitted symbol of s c = 0, probability of successful 
decoding for s p = 1 and s p = are computed as follows: 

Ps alSa , Sp (sc = 0\s c = 0,s p = l) = P(Y^ >t)P(Y^ <t), 

P 8 e \s c ,s p (sc = 0| Sc = 0,s p = 0) = P(Y (Al) > r). 

For the proposed binary MOCSK (BMOCSK), if the 
molecules of type A is used in the current slot, the probability 
of successful decoding can be computed as follows: 

^ c | Sc (sc = 0| Sc = 0) = P(F (A) <t). 

These conditional probabilities for quaternary and higher mod- 
ulation level cases can be calculated similarly. 

III. Performance Analysis 

A. Probability of error 

Let E = l(S c ^ S c ) be the error event. For BCSK, the 
average probability of error is equal to P e = 0.5P(E\S C = 
0) + 0.5P(E\S C = 1). As the decoding of the current symbol 
is dependent on both the last transmitted and decoded symbols 
in BCSK, we have: 

P(E\0)= p s p ,s p \sMp^p\Sc = 0)P(E\S c ^O,s p ,s p ) 
= P Sp( s p) p s p \s p (' S p\ s p) P ( E \ S ^ 0,s p ,s p ), 

Sp . Sp 

P(E\1)= p s p ,s p \sS s P^p\^ = ^) p (E\S c = l,s p ,s p ) 

Sp > S p 

= E p Sp(s P )P §p]Sp (s p \s p )P(E\S c = 1,Wp). (1) 

When s p ^ s p , Pg | S (s p |s p ) is the probability that the 
previous symbol s p is falsely decoded to s" p and can be 
substituted by P{E\s p ). So when s p = s p , Pg , s (s p \s p ) will 
be equal to 1 - P(E\s p ). Therefore, P(E\0) and P(E\1) can 
be computed recursively from ([TJ which is a system of linear 
equations. For QCSK (and higher) levels, the probability of 
error can be computed similarly. 

For the MOSK, the decoding of the current symbol is in- 
dependent of the previous decoded symbol, i.e. Sp. Therefore, 
P e is obtained as 

P e = P(E\s c ,s p )P Sc (sc)P Sp (s P ). 



For the proposed MOCSK, since the decoding of the current 
symbol is independent of both the transmitted and previous 
decoded symbols, P e is computed as follows 

B. Lower bound on probability of error of molecular commu- 
nication system 

In order to determine how well the modulation schemes 
considered above are, we examine a general molecular com- 
munication system with two molecule types A\ and A 2 and a 
memory of only one symbol at the decoder (since we are con- 
sidering the class of simple and practical symbol-by-symbol 
decoders only). Let Bi,l < i < k be the i th input symbol 
(assumed to be i.i.d. and uniform rv's over a set B). The 
encoder uses B\ 1 B 2 , P>i to determine the diffusion rates 

and x[ 2 ^ of molecules of type A\ and A 2 respectively 
at the i th time slot. Thus the actual number of molecules that 
are diffused would follow independent Poisson distributions 

(i) 

with parameters X\ , j = 1,2. Using the thinning property 
of Poisson, the distribution of the numbers of molecules of 
type Aj, j = 1, 2, received at the receiver in the i th time slot 

is - PoissontPiX^ + P 2 X ( i j _\ + A ), where P u P 2 
and Ao are defined in section II-C. The current input symbol is 

decoded as Bi based on the ,Y± and £?;_i (the previous 
decoded symbol). The average probability of error is: 

i k 

p e = -j2 p ( B ^ B *)- 

i=l 

Using the following lemma, a lower bound on P e is derived. 
Lemma 1. I(Bi,Bi) is bounded from above as follows: 

Ifr; %) < I(Y^;XP\X^{) + I(F/ 2) ; X? \xf\). (2) 

Proof: The following chain of inequalities proves d2l. 

I{Bi;Bi) < a I(Bt; B 4 _i, B t - U Y W ,Y} 2) ) 

< I(B i; Bi-i,Bi-i, Y^\Y} 2 \X^\,X^) 
=>I(B ii B i - U Y?\Y?\X&,xV> 1 ) 

= c I{B i -Y?\Y?\B i ^X^X?\) 

<I(xP,X?\B«Yf 1 \YV\B i - U X<%xV) 

</(X« 1 >lf»,B„B,_ 1 ;^,^|xW,X<!») 

= i(xV >x v,yP,yW\ x v >x V} 1 ) 

+ m,B t . 1 -Y^\Y^\X^X^X^\xP) 
= d I(X?\x^^\Y^\X^,X^) 

= H (y?\yV\xI%x(%) 

-H<Y?\y™\xV,X™,X?\x™) 
<H{Y^\x\%xf\) + H[YP\xf\,Xr\) 
- H(Y^\X^,X^) - H(YW\X£\,XP) 

= j(y«;xf)|x«) + 7(y( 2 ) ; xf)|xa). 

Inequality (a) holds because Bi is function of 
-Bj_i, , Y^ , (b) holds as does not provide further 

information about Bi conditioned on X^ x , X^ t , -Bj_i. 
Equality (c) holds because Bi is independent of previous 
message bits, (d) holds because of the Markov chain 
(B u Bi-x) ~ X^X^X^X^ - V«yf>. 
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Using the well-known trick of taking a random index Q 
uniformly distributed on {1.2,- •• ,k} and independent of 
previously defined rv's, one can show P e = P(Bq ^ Bq), 
and the statement of lemma 1 for random index Q instead of 

fix index i. Considering constraints < Xq 



< r r , 



and 



E[X, 



< r av using the upper bounds on the discrete Poisson 



channels derived in [9] if r av /r max £ [1/3, 1] we have: 



;! ^ j)| ^-i)<5 lo g2A7+^log 2 f +Op 17 (1), (3) 



and if r av /r. max £ (0, 1/3) 

< 0.51og 2 PxTmax — (1 ~ 



log 2 (0.5 



o 



Pir„ 



(4) 

)fx — 0.51og 2 2ne 
(1), 

tends to infinity 



fmax 

where Op irmax (l) tends to when Pjr. 
and fx is the solution of r av /r max ■- ^ - ^ erf{ ^y 
Denoting the RHS of the upper bound inequality Q or Q 
by G, and using the upper bound in dZk we have, 



max 
1_ 
2), 



I{B Q ; Bq) = H{B Q ) ~ H(Bq\Bq) < 2G. 



(5) 



Using the Fano's inequality, H(Bq) = log 2 \B\ and the 
convexity of entropy function, we obtain: 



H(P B ) + P e log 2 (|B| - 1) < log 2 \B\ - 2G, 



(6) 



where \B\ is the size of the alphabet of variable Bi. This 
inequality gives a lower bound on the average probability of 
error. Our analysis also extends to M types of molecules in a 
straightforward way; it is omitted due to lack of space. 

IV. Numerical Results 

For numerical evaluations, we have assumed the distance 
between transmitter and receiver and time slot duration are 
d = 16/im and t s = 5.9 sec. respectively and hence the 
probabilities of hitting as Pi = 0.22 and P2 = 0.04 (unless 
explicitly stated otherwise as in the description of two of the 
curves in Fig 1) JS). Also, we have considered Ao = 20. 

Fig. l(Top) shows the plots of the probability of error versus 
the average transmission power (diffusion rate) per bit for 
binary and quaternary CSK, MOSK and proposed MOCSK 
modulations. This figure shows that BMOCSK substantially 
outperforms the BMOSK, with two types of molecules being 
employed in the two schemes. For the quaternary level, 
as expected, MOSK with more complicated transmitter and 
receivers employing four molecule types, yields a better 
performance compared to the proposed quaternary MOCSK 
(QMOCSK) that uses only two molecules types. The same 
results hold for the radio communication as well, where M-ary 
orthogonal modulation substantially outperforms M-ary ampli- 
tude modulation for M > 4, at the cost of higher complexity. 
Also, Fig. l(Top) shows that MOCSK outperforms CSK 
scheme in binary and higher levels. Further this improvement 
becomes more significant as we increase the hitting probability 
P2. To illustrate this, probability of error of CSK has been 
plotted for P2 = 0.1 in binary and quaternary levels in the 
Fig. l(Top). Increasing P 2 results in more interference from 
the previous symbol, thus degrading the CSK performance. 
Generally speaking, as the decoding of the current symbol in 
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Fig. 1. Top: Average probability of error versus the average transmission 
power per bit. Bottom: Probability of error of the prposed QMOCSK and the 
lower bound versus the maximum transmission power. 

the MOCSK is independent of both the previous transmitted 
and decoded symbols, this scheme outperforms the CSK and 
MOSK schemes in a fair comparison. 

Fig. 1 (Bottom) compares the probability of error of the 
proposed BMOCSK and the lower bound derived in ([6]) 
versus the maximum transmission power. In the BMOCSK, 
r av/ r max = 0.5 and 2 molecule types are used. For a fair 
comparison, the lower bound has been computed for M = 2, 
rav/rmax = 0.5 and \B\ — 2. It is observed that the gap 
between the lower bound and the proposed scheme increases 
rapidly. To see how the lower bound changes with increasing 
\B\, the lower bound has been also plotted for \B\ = 4, 8, 16 
and M = 2. 
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